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Modeling Human Interindividual Variability in Metabolism 
and Risk: The Example of 4-Aminobiphenyl 

Frederic Yves Bois, 1 Gail Krowech, 1 and Lauren Zeise 1 


We investigate, through modeling, the impact of interindividual heterogeneity in the metabolism 
of 4-aminobiphenyl (ABP) and in physiological factors on human cancer risk: A physiological 
pharmacokinetic model was used to quantity the time course of the formation of the proximate 
carcinogen, N-hydroxy-4-ABP and the DNA-binding of the active species in the bladder. The 
metabolic and physiologic model parameters were randomly varied, via Monte Carlo simulations, 
to reproduce interindividual variability. The sampling means for most parameters were scaled from 
values developed by Kadlubar et al. (Cancer Res., 51 : 4371, 1991) for dogs; variances were 
obtained primarily from published human data (e.g., measurements of ABP N-oxidation, and ar- 
ylamine N-acetylation in human liver tissue). In 500 simulations, theoretically representing 500 
humans, DNA-adduct levels in the bladder of the most susceptible individuals are ten thousand 
times higher than for the least susceptible, and the 5th and 95th percentiles differ by a factor of 
160. DNA binding for the most susceptible individual (with low urine pH, low N-acetylation and 
high N-oxidation activities) is theoretically one million-fold higher than for the least susceptible 
(with high urine pH, high N-acetylation and low N-oxidation activities). The simulations also 
suggest that the four factors contributing most significantly to interindividual differences in DNA- 
binding of ABP in human bladder are urine pH, ABP N-oxidation, ABP N-acetylation and urination 
frequency. 


KEY WORDS: 4-Aminobiphenyl; toxicokinetics; interindividual variability; population heterogeneity; Monte 
Carlo Simulations. 


INTRODUCTION 

Human heterogeneity is not usually addressed in 
carcinogenesis modeling. Marked interindividual differ¬ 
ences in the activation and detoxification of chemical 
carcinogens 11 ) suggest, however, that certain individuals 
may be at significantly greater risk of chemically in¬ 
duced cancer than the average individual. Increased sus¬ 
ceptibility could be considerable when multiple meta¬ 
bolic pathways and critical interindividual differences in 

1 Reproductive and Cancer Hazard Assessment Section, Office of En¬ 
vironmental Health Hazard Assessment, State of California—Envi¬ 
ronmental Protection Agency, 2151 Berkeley Way, Annex 11, 
Berkeley, California 94704. 


physiologic variables are involved. In this study, we 
modeled the bioactivation, detoxification, and other crit¬ 
ical processes involved in bladder cancer induction by 
the arylamine 4-aminobiphcnyl (ABP), including the in¬ 
dividual variability in formation and activation of the 
proximate carcinogen, N-hydroxy-4-aminobiphenyl (N- 
OH-ABP), and in the fraction of ABP bound to DNA. 

The proposed pathways for human metabolism of 
ABP are schematically described in Figure 1. Activation 
of ABP to N-OH-ABP is mediated by cytochrome P- 
450. (M) N-acetylation, catalyzed by arylamine N-acetyl- 
transferase (NAT), acts as a detoxification pathway for 
ABP induced bladder cancer' 5 ' and leads to the formation 
of N-acetyl-4-aminobiphenyl (N-Ac-ABP). Other detox¬ 
ifying pathways may include cytochrome P-450 medi- 
205 V- 
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Fig. 1 Model of distribution and metabolism of ABP in humans. The 
major pathways for carcinogenesis are indicated by thicker lines. Dot¬ 
ted lines indicate excretion by intermittent voiding of the bladder. All 
processes represented by arrows are modeled by first order linear dif¬ 
ferential equations with rate parameters given in Table I. 


ated ring hydroxylation/ 67 N-glucuronidation (7 ' 8) and 
N-Sulfation/ 97 

The proximate carcinogen, N-OH-ABP, is also me¬ 
tabolized in the liver, forming a N-OH-N-glucuronide 
conjugate 007 and possibly sulfate 01-137 and N-acetyl con¬ 
jugates.^ 047 Some N-OH-ABP binds to hemoglobin (EH?) 
to form an Hb-N-OH-ABP adduct. 05 ' 161 Unbound N-OH- 
ABP is believed to be activated in the bladder under 
mildly acidic conditions via the formation of an electro¬ 
philic intermediate which binds directly to urothelial 
DNA and presumably initiates the carcinogenic pro¬ 
cess/ 4 ' 177 Hydrolysis of the N-OH-N-glucuronide conju¬ 
gate may also contribute to the level of N-OH-ABP in 
the bladder. 

Many of the above activation and detoxification 
processes are mediated by enzymes exhibiting wide var¬ 
iation in activities. Heterogeneity in human cytochrome 
P-450IA2 activity, the isozyme involved in arylamine N- 
oxidation/ 67 has been reported in both in vitro and in 
vivo studies/ 6,18 207 Genetic polymorphism in N-acetyla- 
tion is well established 1 * 1 ' 227 and low acetylation has been 
associated with bladder cancer in a number of stud- 
ies. 15 - 21 - 237 A study by Cartwright, 047 where occupational 
exposure to arylamine carcinogens has been docu¬ 
mented, found nearly all exposed bladder cancer patients 
to be slow acetylators. Theoretically, interindividual var¬ 
iation in sulfation and/or glucuronidation 05 * -277 might also 
contribute to differential hepatic clearance and accu¬ 
mulation of N-OH-ABP in the bladder. 

Finally, two physiological factors, urine pH and fre¬ 
quency of urination, appear to be risk factors for ABP- 
induced bladder cancer. Urine pH may determine the 
rate of formation of the reactive electrophile as well as 


the extent of hydrolysis of the N-glucuronide/ 10) In dogs, 
frequency of urination has recently been suggested to be 
a critical factor in bladder exposure to N-OH-ABP. 12 *) 
Because of the interdependence of these various 
processes, it is difficult to estimate the combined effect 
of their individual heterogeneity on human heterogeneity 
in ABP cancer susceptibility. To describe the possible 
extent of heterogeneity and to examine which factors 
may most contribute to it, we simulated the variability 
among humans by applying Monte Carlo simulations, a 
standard statistical technique/ 311 ' 307 to a pharmacokinetic 
model. Input parameters to the model were derived from 
the literature on human in vitro studies, or from dogs, 
in the absence of human data. ABP carcinogenesis is 
species-dependent and the dog appears to be the best 
model for humans. 1317 The analysis is clearly constrained 
by the limited understanding and data on ABP phar¬ 
macokinetics in humans and on differences within pop¬ 
ulations. 


METHODS 


Pharmacokinetic Model 

Our human model (Figure 1) is similar to that used 
by Kadlubar et al. (2S) to mathematically describe ABP 
metabolism in dogs. All transport and most metabolic 
processes are described by first order kinetics. The body 
is divided into a central compartment and the bladder. 
ABP is absorbed through the gut wall. It then distributes 
to the central compartment (blood and tissues, in Figure 
1) and can be excreted, unchanged, to the bladder com¬ 
partment where efficient reabsorption may take place. 
ABP metabolism, the main elimination pathway, results 
in N-OH-ABP, N-Ac-ABP, or other metabolites, labeled 
ABP-X in the model. N-OH-ABP, in turn, binds to he¬ 
moglobin and other macromolecules, or may be trans¬ 
formed to N-OH-ABP-N-glucuronide, N-Ac-N-OH- 
ABP, or various metabolites (lumped together under the 
label N-OH-ABP-X). N-OH-ABP can also be trans¬ 
ported to the bladder, where reabsorption can occur. A 
similar excretion/reabsorption pathway exists for N-OH- 
ABP-N-glucuronide which can also be hydrolyzed to N- 
OH-ABP, both in the central compartment and in the 
bladder. In urine this process is pH dependent/ 287 Fi¬ 
nally, N-OH-ABP can be activated in the bladder and 
bind to urothelial DNA, a pH-dependent process both in 
rate and extent/ 107 The bladder empties periodically and 
its emptying is modeled as an instantaneous process. 
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Table I. Human Parameter Distributions Used in the Model 


Parameter 

Geometric 

mean" 

CV(%) 

Absorption 

K_ABP_abs c 

0.033 

50 

K—ABP—urine->centraP 

0.033 

20 

K_N-OH__urine—^Central* 

0.011 

20 

Distribution/metaboli sm 

K—ABP—> Ac 8 

0.011 

186 

K_ADP—lOirt 

0.0083 

180 

K_ABP—»X* 

0.0054 

50 

K_N-OH-ABP—>A(V 

0.0027 

186 

K_N-OH-ABP— 

0.0083 

50 

K_N-OH-ABP—>Hb t 

0,0065 

5 

K-N-OH-ABP->Qlu‘ 

0.0036 

50 

K_ hydrolysis—centra]* 

0.0014 

20 

K_ hydrolysis—urine" 

795.0 

5 

K_activation-' 

670.0 

5 

fraction-bound 1 

0.0098 

20 

pH—urine-* 

6 

10 

Excretion 

K_ABP-Hirine* 

0.00039 

20 

K—ABPX—exc* 

o.on 

20 

K _N-OH-ABP—Hirintp 

0.0018 

20 

K_N-OII_N-glu—Hnine-’’ 

0.00036 

20 

Urination Interval 

270 

l 

Time of 1st urination 

135 

m 


* Units; fates in min -1 , times in min. 

* This parameter was scaled from dogs' 28 ) to humans using the allo- 
metric equation X H = X D X{UW H /BW D )-"- , < M1 , with BW„ = 9 kg and 
BW„ = 70 kg. 

' The geometric mean was set at 2/3 of the “ABP—>ABP-M” rate used 
for dogs' 2 * 1 , after scaling to human. At this value about 50% of ABP 
was acetylated. The CV was derived from human data' 3 ’) (sec teat). 

* The geometric mean was set at half of the sum of the K_ABP-»Ac 

and K—ABP-»X rates, corresponding to 33% hydroxylation. The CV 
was derived from human data * 1 ' 1 (see text). 

■ The geometric mean was set at 1/3 the conversion rate of ABP to 
“other** metabolites for dogs' 263 , after scaling to human, correspond¬ 
ing to 17% conversion of ABP to ABP-X metabolites. 

/ The geometric mean was set at one fourth of the “N-OH- 
ABP—»ABP-M” rate for dogs' 28 ), after scaling to human. The CV was 
derived from human data ' 371 (see text). 

* The geometric mean was set at three fourth of the rate of conversion 
of N-OH-ABP to “other” metabolites for dogs' 2 * 1 , after scaling to 
human. 

* K hydrolysis urine (K) was obtained by estimating the rate con¬ 
stant of the hydrolysis (H + and OH - mediated) of N-OH-ABP-N- 
glucuronide (i.e., fitting the equation % hydrolysis = 1 - exp[-K X 
1 q-ph 4 . iq-pok) x time]) to published data' 10 ' * h “ *>. A reaction in¬ 
cluding only H* does not describe the data adequately. 

' K_activation (K) and fraction-bound (F) were obtained by fitting 
the equation /..bound - F X (1 - pH X exp(-K X IQ -' 4 X time) to 
published data' 10 ’ ‘' H ‘ H *1 This equation corresponds to the reactions H‘ 
+ N-OH-ABP -> N-OIT-ADP and N-OH'-ABP + DNA -> DNA- 
N-OH’-ABP, assuming that the second reaction is fast and that only 
a fraction F of the N-OH-ABP is available for reaction. This equation 
for %bound gave the best fit of several alternatives tried. 

> Derived from fitting a lognormal distribution to human data from ,! A 
with truncation interval [4.5, 7.5]' 5J) . 


207 

To account for the human heterogeneity in physi¬ 
ological and biochemical characteristics, the model pa¬ 
rameter values were defined by statistical distributions 
rather than set values. 1 ”'") The specifications for these 
distributions are given in Table I and the derivation of 
these distributions is discussed in the next section. 

Monte Carlo simulations/ 301 were used to obtain the 
distribution of output variables (e.g., fraction of ABP 
bound to DNA in bladder). Each simulation was per¬ 
formed by first picking a random value for each model 
parameter, from its assigned distribution. The oral ab¬ 
sorption of a bolus dose of 100 ng was simulated (our 
conclusions are not affected by the magnitude of the 
dose since linear kinetics are assumed). The time courses 
of output variables, including the area under the quan¬ 
tity-time curve (AUC) of N-OH-ABP and the quantity 
of ABP bound to DNA in the bladder, were recorded up 
to 72 hr. A total of 500 Monte Carlo simulations, the¬ 
oretically representing 500 individuals, were performed. 
Additional simulations were performed to describe DNA 
binding in individuals with average and extreme physi¬ 
ologic and pharmacokinetic status (for Figure 3). 

Characterization of Intcrindividual Variability in 
Parameters 

All parameters, with the exception of the urination 
interval and the time of first urination after ABP expo¬ 
sure, were assumed to be lognormally distributed in the 
population (with mean p and standard deviation SD L ). 
For clarity, Table I gives exp(|l)—the geometric mean— 
and the coefficient of variation (CV) in natural space. 
The SD l used in the computations is related to the CV 
by the approximation exp(SD L )=l+CV/100. For most 
parameters, geometric means were obtained by taking 
the values for dogs given by Kadlubar et al.' 28 - 1 and scal¬ 
ing them to humans. Volumes were assumed to scale 
with body weight and rates with body weight raised to 
the 0.7 power.' 341 The rates reported for dogs sometimes 
apply to a group of metabolites; in such cases the rates 
were apportioned to maintain adequate mass-balance 
(see notes of Table 1). CVs were derived primarily from 
human data, as explained below. For pathways and proc¬ 
esses where no data on ABP or an arylamine congener 


* The geometric mean is the sum of the urinary and fecal excretion 
rates of “ABP-M’*' 281 , after scaling to human. 

' Urination interval was sampled uniformly between the values 60 min 
and 480 min. 

•"Time of first urination was sampled uniformly between 0 and the 
sampled urination time. 
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were available, we used literature data describing the 
human variability in enzyme activities toward a range of 
xenobiotics to estimate CVs. 


Direct Metabolism and Distribution of ABP 

[ % 

!, The variability in the rate of ABP N-oxidation (K. 

' ABP-»OH) was calculated from measurements in hu¬ 

man liver microsomes from 22 organ donors 1 ® and was 
consistent with that of cytochrome P-450IA2 levels in 
hepatic tissue of healthy adults/ 551 Values for ABP N- 
acetyltransferase activities in human liver cytosol have 
been reported in the literature for a small number of 
subjects.™ 1 However, distribution of acetyltransferase 
activities for 2-aminofluorene and the arylamine drug 
sulfamethazine™ 1 are also available and more extensive, 
j The distributions of activities for these two arylamines 

i are very similar in shape and spread/ 35 ’ We therefore 

used the data on 2-aminofluorene activities to obtain an 
estimate of the CV for ABP N-acetyltransferase activity 
(K_ABP—> Ac). 

The minor pathways (ring hydroxylation, N-glu- 
curonidation and N-sulfation) are modeled as a group. 
Based on the data of Butler et al., l6> ring hydroxylation 
j quantitatively contributes very little to the metabolism 

; of ABP. For ABP N-glucuronidation the findings of a 

5-fold range in five individuals’ 81 most likely underesti- 
\ mates the variation due to sm a ll sample size. Other 

findings on the glucuronidation of various drugs indicate 
| wider variation in activities.’ 25 - 27 ' 3 ’™ 25 Data on human N- 

i sulfation of 2-naphthylamine in 20 individuals also sug¬ 

gests substantial variation, with a CV of 40% estimated 
[: from these data. 195 Estimates of the CV for N-sulfation 

of desipramine in humans were even higher. 1435 We 
ji therefore adopted a C.V of 50% for the rate of the com¬ 

bined minor pathways (K_ABP -s- X). 

i 

;j N-OH-ABP Metabolism 

We have designated three pathways for N-OH-ABP 
metabolism: N-glucuronidation, N-acetylation, and X 
which combines all other pathways. Since direct data on 
j N-glucuronidation of N-OH-ABP (rate K__N-OH-ABP 

—» Glu) were limited/’ 05 we used literature data on glu- 
I curonidation of other xenobiotics by humans'' 8 - 27 "' 1 ™ 31 to 

select a CV. There were relatively large variations for 
some xenobiotics (e g., CV > 60%) while for others the 
range in activities was narrow' (CV < 25%). Since the 
specific UDPGT isozyme(s) involved is unknown, we 
assumed moderate variability and selected a CV of 50%. 


For the rate of N-OH-ABP N-acetylation (K_N-OH- 
ABP^Ac) we use the CV calculated for ABP N-acet¬ 
ylation. We assume that pathway X is comprised 
primarily of sulfate conjugation of N-OH-ABP/ 11135 
Given the known genetic polymorphism for phenolsul- 
fotransferase < 2S > we selected a CV of 50% for K N-DH. 
ABP -a- X. 


Other Parameters 

Data from 212 hospital patients’ 101 were used to cal¬ 
culate a CV of 10% for urine pH. The time interval 
between urinations was assumed to be uniformly distrib¬ 
uted between 60 and 480 minutes. The average interval 
(270 minutes) is consistent with published data/ 445 For 
the simulations the time of first urination after exposure 
to ABP was assumed to follow a uniform distribution 
over the time interval between zero and the sampled 
value of the urination frequency. For other pharmaco¬ 
kinetic and physiological parameters, a CV of 20% was 
assumed when variability appeared to be moderate and 
a CV of 5% when little or no variability' was suggested. 
We assumed that FIb-N-OH-ABP binding was mostly a 
physical process and therefore, a CV of 5% was selected 
for the rate parameter K_N-OH-ABP —» Hb. Similarly, 
we assumed that both activation of N-OH-ABP and hy¬ 
drolysis of N-glucuronyl-N-OH-ABP in the bladder 
were purely physico-chemical processes and selected a 
CV of 5% for the corresponding parameters (K_activa- 
tion and K_jhydrolysisL_urine). 


RESULTS 

Figures 2a and 2b are histograms of the AUC for 
N-OH-ABP in the bladder7and of the fraction of ABP 
bound to DNA, respectively, after a simulated oral bolus 
dose of 100 ng. We ran 500 simulations to obtain these 
histograms, each simulation theoretically representing 
the exposure of an individual. The mean and 10th and 
90th percentiles stabilized within that number of runs. 
We ran all simulations up to 72 hr and verified that after 
that time ABP was virtually eliminated from the body 
(in all cases less than 0.5% of the initial dose remained). 
Both distributions are asymmetrical and span 3 or 4 or¬ 
ders of magnitude. The 50% and 95% confidence inter¬ 
vals span approximately 1 and 2 orders of magnitude, 
respectively (the spread is larger for the fraction bound 
than for the AUC). 

A further illustration of the potential differences in 
susceptibility among individuals is given in Figure 3, 
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AUC of the Quantity of N-OH-ABP in the Bladder Fraction of ABP bound to DNA In the Bladder 

Fig. 2 Simulated distribution, for 500 subjects, of (A) the area under the quantity vs. time curve (AUC) for N-OH-ABP in the bladder (units: ng 
X min), and (B) the fraction of initial dose bound to DNA, in the bladder, In (A) the geometric mean is 89.5 ng X min, the arithmetic mean 135 
ng x min, and the 95 pet confidence interval [7.6, 477], In (B) the geometric mean is 0.0034, the arithmetic mean 0.0092, and the 95 pet confidence 
interval [1.3 X 10‘ 4 , 0.053]. 


which presents the time course of the fraction of ABP 
bound to DNA in the bladder, for individuals with high 
and low N-acetylation or N-oxidation activities, and high 
or low urine pH, All other parameters were set to their 
geometric mean. There is a million-fold difference in 
DNA binding between the high pH, high N-acetylation, 
low N-oxidation condition and the other extreme, low 
pH, low N-acetylation, and high N-oxidation. The latter 
extreme differs from the median case by two orders of 
magnitude. 

To rank the various parameters by order of impor¬ 
tance, or their contribution to heterogeneity, the corre¬ 
lations between the Monte Carlo sampled parameter 
values and the model output of interest were exam¬ 
ined.^ 5 ' The correlation coefficients depend on the in¬ 
trinsic influence of each parameter on the outcome 
(which should differ for each outcome), and on the var¬ 
iability of the parameters in the population. The higher 
the correlation, the more dependent the outcome is on 
the parameter. Table D gives the correlation coefficients 
for log-transformed parameters and outputs, and the cor¬ 
responding ranks for the influence of all model para¬ 
meters on the AUC of N-OH-ABP and the fraction of 
ABP bound to DNA, both in the bladder. For the AUC 
of N-OH-ABP the most important parameters are the 

kinetic parameters K_ABP -> OH and K_ABP -» Ac, 

and the urination interval. Figure 4 shows the relation¬ 
ship between the AUC and these parameters. The 
smaller the degree of scatter in the points plotted, the 
greater the correlation and the impact of variability in 
the parameter on the AUC. For the fraction of ABP 
bound to DNA, the three most important parameters are 



Time (min) 


Fig. 3 Simulated time course of the fraction of initial dose bound to 
DNA, in the bladder, for various combinations of parameter values. 

Median conditions:-(K ABP —> Ac = 0.011, K ABP —» OH = 

0.0083, pH 6). Low acetylation, high oxidation, low pH: — + 

(K_ABP -> Ac = 0.00047, K_ABP -> Oil = 0.18, pH 4.5). High 

acetylation, low oxidation, low pH:-x (K—ABP —> Ac = 0.26, K 

ABP —» OH = 0.00038, pH 4,5), High acetylation, high oxidation, 

low pH: — —*— (K_ABP -i Ac = 0.26, K_ABP -» OH = 0.18, 

pH 4.5). Low acetylation, high oxidation, high pH: -+ (K—ABP 

Ac = 0.00047, K_ABP —> OH = 0.18, pH 7.5). High acetylation, 

low oxidation, high pH: v (K ABP —4 Ac = 0.26, K ABP —> OH 

= 0.00038, pH 7.5). High acetylation, high oxidation, high pH: —*— 
(K_ABP -> Ac_f; 0.26, js_ABP —» OH = 0.18, pH 7.5). 


urine pH, and again the kinetic parameters K_ABP —> 

OH, and K_ABP —> Ac (Figure 5); urination frequency 

comes fourth. The highest correlation is between fraction 
bound and urine pH, with a correlation coefficient of 
0.754. 
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Table 11 . Sensitivity" of N-OH-ABP Bladder Exposure and DNA Damage to the Model Parameters 


AUC N-OH-ABP bladder ABP bound to DNA 


Parameter 

Correlation 

coefficient 

Rank 


Correlation 

coefficient 

Rank 

Absorption 

K_ABP_abs 

* -0.0776 

9 


-0.0847 

12 

K_ABP_urine—^central 

-0.0708 

11 


-0.1083 

9 

K_N-OH_urine—^central 

-0.0734 

10 


-0.0416 

18 

Distribufion/mctabclism 

K_ABP—* Ac 

-0.5325 

2 


-0.3808 

3 

K_ABP-iOH 

0-6538 

1 


0.4085 

2 

K_ABP—>X 

-0.1044 

8 


-0.0236 

20 

K_N-OHABP—>Ac 

-0.288 

4 


—0.1961 

5 

K_N-OHABP->X 

-0.1498 

7 


-0.1293 

8 

K_N-OHABP—b 

-0.051 

14 


-0.0258 

19 

K—N-OHABP—iGfu 

0.0545 

13 


0.0862 

11 

K—hydrolysis_central 

0.0208 

18 


-0.0547 

14 

K—hydrolysis__urine 

0.0634 

12 


0.0478 

16 

1 K_activation 

-0.0036 

20 


0 0443 

17 

fraction_bound 

0.0138 

19 


0.1386 

6 

pH_urine 

-•0.0358 

16 


0.7541 

1 

Excretion 

K_ABP->urine 

0.0022 

21 


-0.0090 

21 

, K_ABPX_exc 

0.0302 

17 


0.0514 

15 

K_N-OHABP—nirine 

0.1704 

6 


0.0948 

10 

K_N-OH_N-glu—Hirine 

0.0502 

15 


0.0735 

13 

Time of 1st urination 

0.1795 

5 


0.1314 

7 


“ The sign of the correlations is informative only for the highest correlation coefficients. For the low coefficients the finite sampling (500 runs) can 
lead to values whose signs are not significant (e.g., the coefficient for ABP absorption rate, which should correlate positively with the fraction of 
ABP bound to DNA, is negative). 



Fig. 4 Correlation between the area undeT the quantity vs. time curve for N-OH-ABP in the bladder and (A) the rate of ABP N-oxidation, (B) the 
rate of ABP N-acelylation, (C) the time interval between unnations, for 500 simulated subjects. The correlation coefficients are 0.654, —0.533, 
and 0.364, respectively. 


DISCUSSION 

These results show the importance of considering 
variability from multiple sources when assessing the sus¬ 
ceptibility of individuals to ABP-induced cancer. These 
sources combine in a complex fashion and can be stud¬ 


ied in detail by simulation. As can be seen from the 
results obtained, large differences can arise among in¬ 
dividuals when a number of independent sources of var¬ 
iability exist. These results are consistent with studies of 
cigarette smokers where ABP-Hb adduct levels were 
lowest in individuals who were both fast acetylators and 
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Fig. 5 Correlation between the fraction of initial dose bound to DNA in the bladder and (A) the rate of ABP N-oxidation, (B) the rate of ABP N- 
acetylation, (C) urine pH, for 500 simulated subjects. The correlation coefficients are 0 409, —0.381, 0.754, respectively. 


slow/intermediate oxidizers. 1 "' For ABP, the formation 
of DNA-adducts in the bladder of the most susceptible 
individual is ten thousand times higher than for the least 
susceptible, in a sample of 500 simulations. Whether or 
not this translates into a four orders of magnitude dif¬ 
ference in cancer risk is unknown. However, given the 
probable involvement of DNA lesions in the carcino¬ 
genicity of ABP, it is likely that large differences in 
DNA binding of ABP among individuals correspond to 
large interindividual differences in cancer susceptibility. 

In addition to oxidation or acetylation activities, the 
physiological factors frequency of urination and urine 
pH were observed in these simulations to affect dramat¬ 
ically the extent of DNA binding. Frequency of urination 
appears to be very important in experimental studies in 
dogs as well/ 16 - 28 ! Conversely, some processes were not 
observed to impact the binding of ABP to DNA. For 
example, the reactions involving N-OH-ABP N-glucu- 
ronide were found to be relatively unimportant. This 
could be explained by the fact that, according to our 
simulations, most of the compound is hydrolyzed in 
blood before having a chance to be eliminated in urine 
and therefore is not found in large quantities in the blad¬ 
der. 

Our predictions of interindividual differences, al¬ 
though large in value, may in fact understate the true 
human variability in ABP carcinogenesis since we in¬ 
vestigated only a small number (even if significant) of 
the processes involved. Factors which could contribute 
to either over or underestimation of variabilities are 
model misspecification, unknown cross correlations of 
the parameters, or the mixing of measurement errors and 
intra-individual variability with inter-individual variabil¬ 
ity in the cross-sectional surveys from which the human 
data were obtained. In addition, data are available for 
characterizing variability on relatively few individuals. 


The percentage of ABP metabolized by each path¬ 
way may be a subject of debate, and for many para¬ 
meters where no data are available we took reasonable 
values. Minor pathways of metabolism may be much 
less important than we have assumed. Finally, the ex¬ 
istence of many factors which have not been considered 
here, such as heterogeneity in DNA repair capacity,' 47 - 48 ' 
or the presence of a low acetyltransferase activity in hu¬ 
man urinary bladder cells,' 49 ' may contribute further to 
the variability in the process believed to initiate ABP- 
induced bladder cancer. Had these factors been included, 
the estimated variability would had been higher. 

The present data base to support this type of anal¬ 
ysis for most chemicals, including ABP, is clearly lim¬ 
ited. Human data are available for some critical 
metabolic pathways (e.g., oxidation) for ABP, but the 
numbers of subject studied are too small to give any 
confidence that the ranges reported are representative of 
the general population. The distributions of enzyme ac¬ 
tivities and of most physiological factors were assumed 
to be log-normal. This is a reasonable approximation for 
physiological variables, but the precision of the method 
depends primarily on the quality and amount of data 
available to define the empirical distributions of enzyme 
activities. The sample size of the populations investi¬ 
gated (22 subjects for cytochrome P-450 activities' 6 ' and 
35 subjects for N-acetyltransferase activities' 37 ') is insuf¬ 
ficient to define with precision the shape and the tails or 
extremes of the distributions. The distributions below the 
5th percentile and above the 95th percentile are poorly 
defined and statements regarding the most susceptible 
1/10000 of the population would be inappropriate. The 
number of Monte Carlo simulations performed also in¬ 
fluences the precision of the results, but in this case ad¬ 
ditional Monte Carlo simulations cannot make up for the 
lack of basic information on the distributions of enzyme 
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activities. Extending the simulations beyond 500 runs 
would not significantly improve the determination of the 
extremes or tails of the distributions because the avail¬ 
able data are insufficient for defining the extremes of the 
underlying parameter distributions. Still, the results can 
provide insight on the extent by which risk may vary 
within a population. 

We adapted to the case of humans a pharmacoki¬ 
netic model proposed by Kadlubar et al. r2a) to describe 
ABP pharmacokinetics in dogs. In the new model, acet¬ 
ylation of both ABP and N-OH-ABP, and pH-dependent 
binding of N-OH-ABP to DNA, are explicitly 
incorporated. The model is still a crude approximation 
of reality. For example, inter-species parameter scaling 
had to be used, and the assumption was made that only 
first order reactions were involved. A large literature in¬ 
dicates that cytochrome P-450 and acetyltransferase me¬ 
diated transformations are complex, and may not be well 
modeled by first order reactions/ 22 - 50 ' 52 ' Nonetheless, for 
sufficiently low doses, linearity of the reactions is ex¬ 
pected and the approximation should hold. 

In conclusion, similarly exposed subjects or popu¬ 
lations may have significantly different susceptibility to 
cancer due to differences in the metabolism/activation of 
carcinogens, and differences in physiological status. In 
typical carcinogenesis modeling it has been assumed that 
differences in susceptibility stem primarily from ran¬ 
domness of the occurrence of critical “hits” or DNA 
damage. Our results indicate the importance of consid¬ 
ering inherent differences in human susceptibility to can¬ 
cer as well. Note that most animal experiments, which 
use homogeneous populations provide little basis for as¬ 
sessing human variability. Upper confidence limits on 
risk, when derived from these experiments, are typically 
based on the assumption that each individual faces the 
same risk. Such estimates do not address inherent dif¬ 
ferences in human susceptibility. 
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